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Abstract. 

We present a model to describe the stochastic evolution of stocks that show a strong resistance at 
some level and generalize to this situation the evolution based upon geometric Brownian motion. If 
volatility and drift are related in a certain way we show that our model can be integrated in an exact 
way. The related problem of how to prize general securities that pay dividends at a continuous rate 
and earn a terminal payoff at maturity T is solved via the martingale probability approach. 



1 Introduction 

We consider an ideal model of financial market consisting of two securities: a savings account Z t 
evolving via dZ t = r t Z t dt, where r t is the instantaneous interest rate of the market and is assumed to 
be deterministic (but not necessarily constant); and a "risky" asset whose price at time t: X t , evolves 
according to some stochastic differential eq. (SDE) driven by Brownian motion (BM). As it is well 
known, the prototype model for stocks-price evolution assumes that the return process Rt = \ogX t 
follows a random walk or BM with drift and hence that prices X t evolve via the popular geometric 
Brownian motion (GBM) model, i.e., that X t satisfies 

dX t = nXtdt + aX t dWt ( 1 ) 

Here \i is the mean return rate and a the volatility which are supposed to be constants while Wt is a 
Brownian motion under the empirical or real world probability. We remark that here and elsewhere in 
this article integrals and SDE's are understood in the sense of Ito's calculus. Transition to standard 
(Stratonovitch) calculus can be done if wished. 
The solution to this SDE is given by 

X?™=x exp{vW t +L-^\t} (2) 

After the seminal work of Black and Scholes [1] and Merton [2] , who derive a formula to price options 
on stocks with underlying dynamics based upon GMB, eq. (1) has become the paradigmatic model 
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to describe both price evolution and derivatives pricing. However, while such a simple model captures 
well the basic features of prices it does not quite account for more stylized facts that empirical prices 
show; among them we mention the appearance of " heavy tails" for long values of the relevant density 
probability distributions of returns [3,4]; further, the empirical distribution shows an exponential 
form for moderate values of the returns, which is not quite fitted by the predicted log- normal density 
implied by (2). The existence of self-scaling and long memory effects was first noticed in [5]. Due to all 
this option pricing under this GBM framework can not fully account for the observed market option 
prices and the classical Black-Scholes & Merton (BSM) formula is found to overprice (respectively 
underpricc) "in (respectively, out of) the money options". Apparently, for empirical prices of call 
options to fit this formula an extra dependence in the strike price, the volatility smile, must be 
introduced by hand. 

After the seminal paper by Mantegna and Stanley who studied the empirical evolution of the stock 
index S&P500 of the American Stock Exchange during a five year period, several authors have elabo- 
rated on the possibility that price dynamics involves Levy process and have discussed option pricing 
in such a framework (See [5-13]). For complete accounts of option pricing and stochastic calculus from 
the economist and, respectively, physicist, points of view see [14-17] and [18-21]. 

Here we shall focus in another different aspect that some traded stocks seem to present, viz the 
possibility of having, at some level, strong resistances both from above or below. For example, 
corporations or major institutions may have laid out a policy under which heavy buy orders are 
triggered whenever the stock price hits this level. Such feature can not be described with Eq. (1) as 
under such an evolution prices can reach any value in (0,oo). Concretely, in this paper we want to 
model the evolution of a market stock which has a strong lower resistance at some level c where we 
suppose that c is a constant. 

In section (2) we present a model that incorporates an attainable barrier at the point x = c > and 
hence can, in principle, be used to account for such a fact. We next derive the evolution of the asset 
and the probability distribution function. It turns out that c is a regular barrier in terms of Feller's 
boundary theory and hence a prescription on how to proceed once reached must be given. In section 
(3) we study pricing of securities under such a model and obtain a closed formula for valuation of 
European derivatives that have, in addition, a continuous stream of payments. We tackle this problem 
using the Martingale formalism of Harrison et al [22] and obtain the partial differential equation (PDE) 
that the price of a security satisfies. Solving this eq. corresponding to particular final conditions we 
obtain the price of options under this model. This price is compared with that given by the standard 
Black& Scholes- Merton formula. In the appendix we consider some technical issues concerning value 
of the market price of risk and the the existence of the martingale measure or risk free probability 
under which securities are priced. 

2 Price evolution under the martingale probability 

Let r t be the deterministic interest rate at time t and Z t = exp J Q r s ds be a "savings account " . As we 
pointed out we consider that X t is the i-price of a tradeable asset that has a strong lower resistance 
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at some constant level c where < c < xq = Xq. Mathematically this implies that the values of 
X t must be restricted to the interval [c, oo) and hence X t must have a boundary point of a certain 
kind at x = c. From intuitive financial arguments the boundary can not be of absorbing type since in 
that case, once reached, the price X t remains there. Further it seems reasonable to assume that there 
exists positive probability to attain the boundary; we suppose that this event "triggers" bid orders 
and hence that Xt ricochets upon hitting the boundary. Therefore in such situation the assumption 
that prices evolve via Eq. (1) is no longer valid. The obvious modification wherein prices evolve as 

X ( GBM = c+(x - c )e CTm+ (^-^) t (3) 

is also ruled out as this evolution implies that c < X t < oo but the value X t = c is never attained 
and the probability to get arbitrarily close to the barrier tends to zero with the distance to it. Thus 
trajectories corresponding to a such model never quite seem to reach the support (In terms of Feller's 
theory briefly reminded below x — c is a natural barrier at which Feller functions blow up). 

Motivated by similar ideas in the context of the Cox-Ingersoll-Ross model of interest rate dynamics 
[23] we now introduce a more satisfactory model which satisfies the aforementioned features and is at 
the same time analytically tractable; we shall suppose that X t evolves via the SDE 



dX t = a(t, X t )dt + b(t, X t )dWt, where a(t, x) = fix, b(t, x) — a(t)\/ x 2 — c 2 



(4) 



Here Xq = xq > c, /i is the stock mean rate of return and b(t,x) the volatility coefficient. Indeed, 
under such a dynamics it follows from (4) that as x approaches the point c, b(t, x) tends to zero and 
hence X t evolves roughly like dX t = /iX t dt implying that X t will increase and then escape from the 
boundary. 

For valuation purposes one needs to consider the evolution under a new probability that might be 
different to the empirical observed probability. Mathematically speaking a such a probability is defined 
requiring that under it the discounted prices X' t = Xt jZt are martingales (this risk-neutral probability 
was introduced in [22] although the underlying idea pervades the original work of Black-Scholes & 
Merton [1,2]). Stated another way, this means that under the risk-neutral probability, the stock price 
X t evolves, on average, as the riskless security Z t thereby preventing arbitrage opportunities. Indeed, 
the martingale property implies 



Xn) = E* 



x o) = x o = xo 



(5) 



where E* {X t 
Hence 



Xq ) is the conditional average of X t given Xq with respect to the martingale probability. 



E* X 



X I = x Z t 



(6) 



More generally, given the past history !F S of the process up to time s (i.e., the cr-ficld of past events) 
one has 



E* X 



Z f E* XI 



F s ) = ZtX' = X s eJ>' 



(7) 



We shall assume that our market is efficient, i.e., that the martingale probability P* exists- which is 
not always the case. In such a case the explicit form of the original drift coefficient is only needed 
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to go back to the empirical or real world probabilities. Indeed, it follows from these arguments that 
consideration of this probability amounts to redefining the evolution equation without changing the 
volatility coefficient b(t,x) but replacing the drift coefficient to a*(t,x) = r t x, independent of the 
initial coefficient a(t, x) = /jlx. 

Unfortunately, in general it is not possible to solve the redefined SDE corresponding to the diffusion 
coefficients (4) with a* (t, x) = rtx. However, it turns out that in the particular case when a 2 (t) = 2r(t), 
i.e., for eq. (8) below, then both this SDE and the prizing problem can be solved as we next show. We 
shall consider this case and hence we suppose that under the risk neutral probability P*, X t evolves 
via the SDE 

dX t = nXtdt + ^2r t (X 2 ~ c 2 )dW t *, X = x >c (8) 

Here Wj* is a BM with respect to the risk neutral probability. Technically, the existence and nature 
of all objects introduced below is a difficult point. In the appendix we sketch how to perform such a 
construction. 

In the sequel all quantities are referred to the probability P* and hence X t evolves via (8) which is 
our fundamental equation. Further for ease of notation we drop here and elsewhere the use of *. 
The return process Rt = \ogX t /xo is found via Ito's rule to satisfy 

dR t = r t c 2 e- 2Rt dt + ^2r t (l ~ c 2 e- 2R ^jdW u R = (9) 
Thus only when R t is close to it does behave like a classical random walk. 

Useful information about the behavior of the price process at x = c follows by careful inspection of 
the nature of the boundary x — c. Consider the Feller functions S(c, x), f2(c, x) defined by 

where p(x) = \J x 1 — c 2 . The reader is referred for these matters to [24]. Notice that the integrand is 
singular since it has a square root singularity at x — c. Upon evaluation of the integrals we find that 

S(c, a) = Q(c, x) = 1 log 2 ( £±^Eg ) (ii) 
4r c 

Thus, unlike what happens in the model defined by (3), we have that S(c, a;) = f2(c, x) < oo are finite, 
corresponding to a regular boundary which can be both reached and exited from in finite time with 
positive probability. 

While Feller analysis shows that the boundary is attainable it does not clarify if the process can 
be continued past the boundary (and hence whether prices below the level x = c can be attained). 
Further it it is unclear what is the probability to reach the boundary or how the to continue the 
process upon hitting the boundary. 

These kind of problems regarding behavior of the process at and past the boundary are generically 
quite difficult to tackle. It turns out that for our particular model the behavior of the process is 
completely determined. Actually we have found that the solution to eq. (8) is given in a fully explicit 
way by 
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X t = c cosh ( J V2r~ s dW s +k) (12) 

where k = cosh - (xo/c). A typical path t — > X t is shown below in Fig. 1. To prove (12) let X t = g(Y t ) 
where 

g(z) = ccosh (z), Y t = n+ f s/2r~ s dW s (13) 

Jo 

Using Ito's rule and the fact that dY t = \J2r t dWt we find that X t has diffusion coefficients a(t, x), b(t, x) 
satisfying at x = g{z) 

a{t,x) = ^-d^g = r t g(z) = r t x 



b(t, x) = V2F t d z g = V2F t c sinh (z) = J 2r t [x 2 - c 2 j (14) 
i.e. X t = g(Y t ) solves the SDE (8). 

Notice that the last equality and the fact that the sinh takes both signs, imply that the following 
prescription must be given at the barrier: 

b(g(z)) = -b(g(-z)) (15) 

Thus (12) solves (8) provided the square root is defined with a branch cut on (c, oo). 
We note that the conditional density f(s, Y\t, y),t < s of Y t solves 

(d t + r(t)d yy ) / = 0; f(s, y\s, y ) = 5(y - y ) (16) 

which is converted into the classical BM or heat equation upon time transformation via t' = ip(t) 
where we define 

<p(q, r) = 2 / r s ds; <p(t) = <p(0, t) = 2 / r s ds (17) 

It follows that the process Y t has the distribution of a BM evaluated at time t' and hence we can 
represent X t as 

X t = ccoshF t where Y t = n + B^u) (18) 

and Bt is a new BM. It follows from (12) that Xt > c and that Xt attains the barrier c whenever the 
process Y t reaches , i.e., when J* \/2r s dW s + n = 0, which, according to (18) happens eventually 
with probability one. As pointed out, it follows from (18) that in that case the process X t is reflected 
and hence the level x = c > acts as a resistance of the stock value. 

Let us now obtain p(T, X\t, x), the probability density function (pdf) of the price process conditional 
on the value at time t: Xt = x,t < T. This pdf satisfies the backwards Kolmogorov-Fokker-Planck 
equation 

(d t + r(t)(x 2 - c 2 )d xx + rxd^p - 0; p(T, X\T, x) = S(x - X) (19) 

Motivated by (13,17) above we define new coordinates t' — ip(t), y = log (x + \/x 2 — c 2 ^j — logc. In 
terms of the new coordinates p solves 

(d t + ^d vy y = 0;p(T,Y\T,y) = S{ccoshy-ccoshY) (20) 
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Using the well known formula 

cS(c cosh y — ccoshF) = (%-Y) + <% + y))/|sinhy| 

we find that 



p(T,X\t, x) = - 1 ]T exp ( - — * log 2 " + ^IZ_ ) (21) 



We next compare the evolution of prices under this model and that described by GBM. For a mean- 
ingful comparison we need to have r t — r constant and, in Eq. (1), fj, = r,a = s/2r. In this case (2) 
yields 

X? BM = x Q e^~ rWt (22) 
Note that whenever Wt >> 1 both process behave in a very similar way: 



lim X t /X? BM 



x + ^fxl - 6 
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Wt^oo 2xo 
However as Wt — > — oo then 

X t « (x - yJx 2 Q -<?)/2e^ r \ w '\ -» oo, X? BM = x e^ rW * -» 
Finally, we note that in the limit c — > we can recover (2) from (18). A careful calculation shows that 

limX t = xq exp £>2rt 
and hence, using that y/2rWt is a BM at time 2rt, we recover (2). 

3 valuation of securities 

We consider here the valuation of securities earning a terminal payoff Qi(Xt) at maturity T. We also 
allow for the security to pay dividends at a continuous rate 02(s, X s ),t < s < T where we suppose 
that both 0i : K — * K and 02 '■ R + x 1 ^ I are continuous. The standard case of stock option 
valuation corresponds to taking Qi(x) = (x — k) + , 02 (x) = where fc is the strike. 
Let ft be the (actual) i-price of such European derivative maturing at T, which must also depend on 
T and the actual price x = X t of the stock. 

We assume the existence of risk-neutral probability P* under which relative prices of stocks and more 
generally, of self-financing strategies v' t = Vt/Z t are martingales with respect to the history of the 
process up to time t: Tt (notice that we shall drop again the symbol*). If this is the case, reasoning 
similarly as in (7) and use of the martingale property yields that 



v t e^ ridl = Z T v' t = Z t e(v' t T^j = ¥.(v T 
To continue further we note that vt must satisfy at t = T that 



v T = &i(X T ) + 



j e 2 (s,X s )Z T /Z s ds = v^ + VT (25) 
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as the RHS is precisely the earning at maturity. Note further that = Oi (Xt) is obviously Markovian 
and hence that E^vt Ttj = E^«t Xtj ■ 

It can be proven that v T = J. T & 2 (s, X s )Zt/ 'Z s ds is also a Markov process and hence that it satisfies 



E(^ e 2 {s,X s )Z T /Z s ds^ t ) =E(J <d 2 {s,X s )Z T /Z s ds\x t ) 
Hence we finally obtain the price of the security as 



Vi 



e-f* Tr - ds G 1 (X T ) + £ e 2 {s,X s )Z t /Z s ds 



X, 



(26) 



(27) 



This can be simplified further by reasoning as follows. Let Xl> x ,t < s be the price process at time s 
knowing that it starts at x at time t, i.e., Xt' x = X s is the solution to (8): 

dX s = r s X s ds + \/2r s [X 2 — c 2 ) dW s , with initial condition X t = x. If we use the well known 
property ([24]) 

X t = x] = E(x^ x ) (28) 



E(X T 

then eqs. (13,27) are rewritten in the convenient form 

X l T ,x = ccosh (BpfaT) + v )i v = cosh~ 1 (x/c), 



E 6i(X T ) 



X t = x) =E [Q 1 (X, 



dY 



9i(ccosh(y + v)) V- 



(29) 



and so forth (actually, X T ,Xt = X^' x ° when r is constant). 
Alternatively, with (21) at our disposal we also obtain that the price of a security paying dividends 
at a continuous rate & 2 {s,r s ) and a terminal value Oi(ry) at maturity is given in an explicit way by 



e -f? ridl J dXQ 1 (X)p(T,X\t,x) + J ds J dXQ 2 (s,X)p(s,X\t,x)e-f* 



ndi 



(30) 



Note that by using the Feynman-Kac theorem vt may be also evaluated by solving the backwards 
equation 

(d t + r(t){x 2 ~ c 2 )d xx + rxd x -r)v= -Q 2 (t, x) 
with the terminal condition limn* = Qi(x)- 

As a natural application we evaluate the price of the plain vanilla call with strike k corresponding to 
Qi (x) = (x — fc) + , 82(2:) = 0. Let us introduce 



x=(x+ \/x 2 -c 2 )/2, k= {k + ^fc 2 -c 2 )/2 
^<p{t,T) k 



N = 



M 



1 



log— ,M ± = M±y/<p(t,T) 



(31) 
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Then, in terms of <I>, the distribution function of the normal variable 7V(0, 1), we find from (27) to 
(30), that, if Xt = x > c, the plain vanilla call price is given by 

Vt = x($(N + ) + $(M_)) + (x - x)($(N-) + ®(M + )) - ke~ J? r W da ($(N) + $(M)) (32) 

The situation when c = and r is constant amounts to having no barrier and hence (32) must reduce 
to the BSM formula. Indeed, in this case one has x = x, (p(t, T) — 2r(T — t), 

N= N° = , 1 logy, Nl = N°± J2r(T -t),M= -oo (33) 
v/2r(T - 1) fc ± V V ' 

and hence $(M) = $(M±) = 0, most of the terms in (32) drop out and we recover the BSM formula 

vt = x$(Nl) - ke-^-'^iN ) (34) 

Another interesting simplification appears when the strike k coincides with c: k = c. Using that for 
this case is k = c, M = —N, M± = —N T and the well known property $(z) + z) = 1 we find that 
all parenthesis in (32) add to 1 and the price of the option is the deterministic price: 

v t = x - ce~ St 'K s ) ds . 

The result is easy to understand; indeed, as we pointed out (12) implies that X t > c which rules the 
possibility to have k < c, i.e., this case corresponds to having taken k at its lowest possible value. 
Thus all uncertainty disappears since with probability one Xt > k and the option will be exercised 
will probability one. 

In figure (1) we plot a typical path of the price process (8) starting at xq = 5. We assume a yearly 
interest rate r = 4.5% yr , annual volatility a = 30% and suppose that the support is placed at 
c — 4. Notice how eventually prices get near and eventually hit the support level c — 4 lingering 
around for some time. 
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Figure 1: A possible evolution of the price process. We plot X t as a function of time during a time 
span of five years (t = 5). The parameters have been chosen as: r = 4.5% yr _1 , a = 30%. 

In figure (2) we plot the call price v in terms of the initial stock price x corresponding to a constant 
annual interest rate r = 4.5% yr " 1 with annual volatility a = 30% and time to maturity T — t = 1 
year. The barrier is located at c = 4 while the strike k = 8. The thick solid line represents (32) while 
the dotted line is the classical BSM call price; the deterministic price d=(x- ke-^ r ^ ds )+ is the 
thin straight line. In all cases one finds that the prices implied by the classical BSM valuation formula 
and (32) are quite similar, specially for long x. Notice how the BSM formula always overprices the 
call option compared with the formula (32). However, it seems that the variation is only significant 
in the region c < x < 2c, irrespective of how large c is. 

Actually, we find in all cases that the deviation of (32) from (34) is quite small (see figs. 2 and 3). 
This is easy to understand qualitatively when x is long, since then x — x — |- + 0(l/x 3 / 2 ) and hence 
we find the expansion 

VW7r)N = io g (|) + bgc - ^ + o(i/* 3/2 ), 

where < C = logfc/fc < log 2. Further M « (1/ logx)e" log2 x and 

v t ^ x - ke~ ^ r( - s)ds + 0{l/x) 

as in the BSM formula. However, such a rough argument does not account for the close similarity 
even for small and moderate values of x. 

The dependence of the option price upon the time to maturity up to 40 years is shown in figure 3. 
The plot corresponds to an ATM option for which moneyness x/k = 1. The rest of parameters have 
not been changed. 
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Figure 2: plot of the call price v in terms of x corresponding to an annual interest rate r = 4.5%, 
annual volatility a — 30% and T — t = lyear with a barrier at c = 4. Strike is taken as k = 8. The 
thick line represents (32) while the dotted line is the BSM call price. Deterministic price is the thin 
line 
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Figure 3: plot of an ATM call option price v in terms of time to maturity z = T — t. The thick solid 
lines represents (32) while the dotted line is the BSM call price. Deterministic price is the thin line. 
Parameters are as in the previous caption. 
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4 Appendix 



Consider again eq. (4) dX t = /iX t dt + b(t, X t )dWt where Wt is a BM with respect to the real world 
probability and b(t, x) = yj2r t (X? — c 2 ). Let 

u{t,x) ee (r(t) -fj)x/b(t,x), M t = exp(-^ u(s, X s )dW s - ^ J u 2 {s,X s )ds} 

Here u(t, x) is the so called market price of risk . Then if X t solves the above equation M t can 
be proven formally to be a martingale. Note however that a rigorous proof of the latter fact runs 
into technical difficulties due to the singularity of u at x = c which might prevent, in principle, for 
M t to be a Martingale. We skip a rigorous analysis as we expect this to be the case. Defining the 
risk neutral probability P* by dP* = MycZP it follows from Girsanov's theorem (see [17,20]) that 
Wf ee Wt + J Q u(s, X s )ds is a BM under P*. In this case an easy calculation shows that X t also 
satisfies Eq. (8) driven by W t *, a BM with respect to P*. 
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